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" ■ I. INTRODUCTION 

q ; 

The solutions of periodic systems with periodic boundary conditions (PBC) and of finite periodic systems with 
complete quantum confinement (CQC), are characterized by a discrete spectrum of the Bloch quasi- momentum q. In 
contrast in finite systems with open boundary conditions (OBC), solutions do exist for any real value of q, as in an 
''O ■ infinite periodic system. 

In the literature general orthogonality relations for Bloch waves can be found only for either (i) infinite systems or 
(ii) systems with either PBC or CQC, but not for systems with OBC. Nevertheless, locally periodic regions of finite 
length L with open boundary conditions are more realistic models of crystals than systems with periodic boundary 
conditions which have been studied extensively in the literature. 

In this paper we derive vanishing integral relations valid for particularly chosen Bloch functions occuring in finite 
periodic systems consisting of an integer number of unit cells with either PBC, CQC or OBC. The latter pertain 
to finite-size crystal models as well as semiconductor superlattices. These integral relations include orthogonality 
relations and conditions for the vanishing of matrix elements. The choice of Bloch functions in these relations is 
directly related to the symmetry group of a system in which periodic boundary conditions are imposed at the ends 
of the finite periodic potential. However, the relations we derive are not restricted to "quantized" values of q which 
r*"** \ correspond to systems with PBC or CQC only, but hold also for uncountably many well-defined combinations of 
wave numbers q from the continuous quasi momentum space. On this basis we study properties of specific linear 
combinations of Bloch functions. Our results indicate that the large body of group theoretical results on systems with 
PBC can also be applied to more realistic models with no apparent symmetry whatsoever. 

The paper is organized as follows. In Section |TT] we shortly review some important properties of the solutions 
of the Schrodinger equation in the case of (finite) periodic one-dimensional potentials. In Section IIIII we state the 
fundamental result on vanishing integral relations for finite periodic systems [cf. Eq. (fT5|) . (|14p and (jTSJ)], which is a 
generalization of the well known selection rule theorem for PBC. Subsequently in Section [IV I we discuss applications 
O ■ of this main result to vanishing integrals involving resonant Bloch functions [as defined in Eq. (jlip j. vanishing matrix 
elements of specific Bloch functions and orthogonality relations. Finally in Section [V] we consider special properties 
- among them time-independence and real valuedness - of physically relevant expectation values of Bloch functions 
and superpositions thereof in open finite systems. Appendix |A1 contains concepts which are subsequently used for a 
group theoretical proof of Eq. (fl"3|) . 
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II. SCHRODINGER EQUATION FOR (FINITE) PERIODIC ONE-DIMENSIONAL POTENTIALS — 

FUNDAMENTALS 



It is well known that the fundamental solutions of the time independent one-dimensional Schrodinger equation with 
an infinite periodic potential (with period d) 



*(£C) = (x) , Vper (X + d) = Vper (x) , 



inside an allowed band are Bloch functions pj ^f(x), given by 



"9? (x) — u q (x) exp(iqx), Uq(x + d) — Uq(x), q G R, 



(1) 



(2) 



where g denotes the energy dependent Bloch wave number or quasi momentum 0] ■ To obtain the Bloch wave number 
q in the "reduced zone scheme" Q the minimal residue of q mod 2ir/d is taken. As is well known this "folds" the 
Bloch wave number from the n-th to the first Brillouin zone q £ (—ir/d, ir/d]. In this reduced zone scheme the original 
Brillouin zone number is added as band index n. Then Eq. ([2]) reads 



(x) = ^n,q( x ) = U n ,q{x) exp(iqx), U n<g (x + d) 



;0)> 



(3) 



which satisfies H^^ q (x) — E n ,q^n q (x). The relation E n (q) := E n . q is known as a dispersion relation. Due to 
Kramers degeneracy E n _ q = E Ut - q . Therefore any superposition [3[ 



with arbitrary satisfies H^f^ q _ q ^(x) = E n ^®^ q) {x). 



(4) 



A. Periodic boundary conditions 



Since the introduction of Bloch functions their properties have been studied in detail mostly for systems with 
periodic (Born- von Karman= BvK) boundary conditions, i.e. vE'(x) = ^{x + Nd). The BvK boundary conditions 
imply a quantization of the Bloch wave number in Eq. ([3]) in the form of 

&=2b/W, kGK BvK :={-[(N-l)/2\,..., -1,0,1,..., [N/2\}cZ, (5) 
but no restriction on a + and oT in Eq. [[3]), The set 

ObvK := {q^ K , k G K Bv k} , |0 Bv k| = N (6) 

forms a finite cyclic group (QbvK,©), where the group operation © is addition modulo 2-7r/d with minimal residue 
taken, i.e. 

{<?i + <?2 + 2n/d if qi+q2<-n/d, 
qi + <72 if - n/d < q 1 + q 2 < ir/d, (7) 

qi+q2-2n/d if qi + q 2 > n j ' d. 



B. Complete quantum confinement in finite periodic systems 

Recently, interest in studying properties of finite locally periodic quantum systems has turned up [3, [H, B, H, El H, 
E3, EH- One option to address such finite systems is to introduce complete quantum confinement (CQC)[(| EI, i.e. 
^/(0) = ty(Nd) = 0. The wave number quantization of the real Bloch wave number in Eq. ^ in this case yields a set 

£cqc := {$b = kir/Nd, k G {-N + 1, . . . , -2, -1, 1,2, . . . ,JV - 1}} , (8) 
Iflcqcl = 27V -2, 

but with the additional restriction: a + +a~ =0, if u nyq (0) is choosen to be real [Hj]. Note that due to this constraint 

the wave numbers <7cqc an< ^ Qcqc arwavs jointly appear in Eq. ((4]). The set Ocqc coincides with the set rcs below 
[cf. Eq. pTjl ] and can analogously be closed to a group. 
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C. Open boundary conditions for finite periodic systems 

The potential of an open, finite TV-fold periodic systems d, 0, H, E3] is 

( ifx<0, 

V(x) = \ V per {x) if < x < L = Nd, (9) 
[ if x > L. 

The (rigged) Hilbert space of this problem is not restricted to the interval [0, L) but extended to the whole real line. 
In an open system, solutions (also called scattering states) in the form of Eq. ([4]), exist for < x < L for all Bloch 
wave numbers in (— n/d, ir/d]: 

!Al exp[ik(E n . q )x] + Bl exp[— ik{E n:q )x] if x < 0, 
a + $B q (x) + a-y%_ q (x) if < x < L = Nd, (10) 

A R exp[ik(E n , q )x] + B R exp[-ik(E n ^ q )x] if x > L. 

Hence the Bloch wave number q is not quantized, and there are no restrictions on the coefficients a^, as in the PBC 
system. 

By considering resonant states (with a transmission probability of unity and a reflection probability of zero) of an 
open system with a potential given by Eq. (J9j) a new kind of "quantization" arises. The set of discrete g-values that 
belong to Fabry-Perot resonances [Hf of the open system is given by 0, 0, H, EH 

Ores := {q^l = kir/Nd = 2kir/2Nd, k G {-N + 1, . . . , -2, -1, 1,2, . . ., N - 1}\ , (11) 
|0 rcs | = 2N-2. 

This quantization corresponds to a BvK system with 2N periods. The "doubling" of the number of periods is 
fundamentally related to the fact that resonant wave functions fulfill the relations [HQ ^(0) = ±^(Nd) which means 
that \&(0) = x S(2Nd). It is possible that additional non- Fabry-Perot resonances exist which are due to transmission 
resonances in the unit cell 0, H, 0, EH but these are outside the scope of the current paper. 

Note that the union of Q TCS with the set of q values belonging to the band edges, i.e. res := Ores U {0, ir/d}, forms 
together with the operation © a group and |0 res | = 2N. Obviously, as 2qiei = q^ K , res = HbvK u Qrll © ObvK 



III. VANISHING INTEGRALS INVOLVING PRODUCTS OF BLOCH FUNCTIONS 

In this paragraph we state fundamental integral identities which emerge in the analysis of Bloch waves in finite 
periodic systems (see e.g. Ref. To our best knowledge these identities are presented below in a general form for 

the first time. 

The group (ObvK,©) [cf. Eq. ^] is a discrete (normal) subgroup of the compact continuous one-dimensional 
rotation group (Qf,©), where 

Q*:=(-7r/d,7r/rf|. (12) 
The following fundamental integral identity holds: 

^3 q r G 0bvk\{O} / dx ^ [ I = 0, (13) 

r— 1 r—1 

where q r G £lf, (J) denotes a sum with respect to the group operation in and f per {x + d) = f per (x) denotes an 
arbitrary d-periodic function. 

Equation (fTB"]) can trivially be extended to 

l§ 9 7 l§ 7 U} X ' ; '' — ^ — 1 1 — ^ — 11 — ^ — =0 ' (14) 

where we make use of the fact that ^„*q an( i -q differ at most by a constant phase factor. 
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Most generally one can extend Eq. (fT"3| to 



[■Nd 

I dxl(x) = 0, 
Jo 



(15) 



where the integrand I(x) is any sequence of terms belonging to the set 

{d/dx,*l >qr (x),{*l tq> (x)]*,fr(x)}, (16) 

where additionally pairs of parentheses can be inserted appropriately and operators act from right to left as usual. 
Once again 

® (® (_<Zs) ) G QbvkU ° } - (17) 

For example valid integrands are dj 'dx^ B iqi (x)ff er '(x)^ B 2m (x)[d/dx^ B 3iq3 (x)]f^ er (x) (with qi®q 2 (Sq 3 e O B vk\{0}), 
where the first d/dx acts on the rest of the term, or qi (x)]*d 2 /dx 2 ff er (x)^ B 2 q2 (x) (with -qi © q 2 G Obvk\{0}) 
where d 2 /dx 2 acts on the product ff er (x)^ B ^ q2 ( x )- This last generalization follows from Eq. (fT4"|) by applying the 
product rule of differentiation. 

The proof for Eq. (TT3|) can be seen from two perspectives. On the one hand, it is well known [2| that 

dxf pcT {x) exp (iq^ K x\ = °> if fc ^ 0, (18) 
and obviously 

J dxf pcT (x)exp ^i^2q r x^ =0, if ^ 1r S ObvK ^ 0, (19) 

for otherwise arbitrary wave numbers q r from £JJ?. Equations (|13p and (|19[) are equivalent taking into account the 
translational invariance of the differential operator. Normally relations such as Eq. (|19p are widely used in solid state 
physics but all summands in q r are restricted to belong to ObvK as only PBCs are considered. Eq. (Tl3|) extends 
the well known selection rules theorem to uncountably many well-defined combinations of wave numbers q from the 
domain of continuous quasi momenta for finite locally periodic open systems. 

In Appendix A we discuss the proof of Eqs. (|13H15|) from a rigorous group theoretic point of view. 

IV. APPLICATIONS TO SPECIFIC CASES 
A. Vanishing integrals involving two resonant Bloch functions 

Identities involving two resonant Bloch functions ^ B (k) and ^> B , kn which follow from Eq. (fT3l) are 

r Nd dU ( ^d l ^ B , k) (x)d m ^ B , h ,Ax) 
Qid © ^ cs ; G QbvK \ {0) dx — — o — t — 0, (20) 

where j, I, m £ IN U {0}. 

For example an extension of Eq. (A2) in Ref. [l(| to position dependent effective masses m{x) is of the form of 
Eq. fH with k = k',j = 0, f per (x) = l/m(x),l = 0,m= 1. 

B. Vanishing matrix elements for specific Bloch functions 

Equation (TT4|) allows to consider matrix elements of the form {^ B t q .\A\^ B . qj )Nd, where A is a periodic differential 
operator with respect to space, i.e. 

d l 

A(x) = ^2ai(x)— -, ai(x + d) = ai(x) (21) 



dx 1 
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and (-)Nd denotes projection onto the finite space region [0, Nd], i.e. 

{V\A\<I>) Nd := (*|AP FPS |$), (22) 

where Pfps = J dx \x)(x\. Since the rotation group [cf. Eq. (JT2J)] can be trivially decomposed in cosets with 
respect to ObvK : 

&c = (J © Hbvk) . Vc Jd € (-Tr/iVrf, Tr/iVd] = (-gW , g&>], (23) 

the following relation holds for any operator A of the form given in Eq. ([?!]) : All matrix elements of A for two different 
Bloch waves (of the same or of different bands) with wave numbers qi and qj in the same coset of 0bvK vanish, i.e. 

Qi + qj\Qi,Qj e q? d © 0BvK => = 0. (24) 

Of course, the condition is equivalent to qj — = 7^ 0. 

A special case of this situation is the vanishing of matrix elements of (Kramers-) degenerated resonant solutions 
which emerge in the fundamental solution of the open or bounded finite periodic system, 

<f B i±u) \A\¥ B (Tk) ) Nd = 0, (25) 

f-fc) (fe) _ (k) .-^ 
Since qi-es ,qres £ qres © UfivK- 



C. Orthogonality relations 

Setting A = 1 in Eq. (|24|) we obtain the following orthogonality relation: Two different Bloch waves (of the same 
or of different bands) with wave numbers qi and qj in the same coset of ObvK are orthogonal, i.e. 

Qi + '/..:</-'/.. G <zf rf ©fl B vK =► <*£, 9 J*£, w }jvd = 0. (26) 

Again, the condition is equivalent to qj — qi = 41? 7^ 0. 

Similarly to Eq. (|25|) we get orthogonality of (Kramers-) degenerated resonant solutions, i.e. 

<* s (±fc) |* B (Tfe) )jv<i = 0. (27) 



V. EXPECTATION VALUES FOR SPECIFIC SUPERPOSITIONS OF BLOCH FUNCTIONS 



We consider (normalized) superpositions composed of Bloch waves from arbitrary bands (with band indices ni) for 
which all reduced Bloch wave numbers qi belong to the same coset of 0BvK but are different, i.e. 

|*> = X) a *l*^,«>' HH 2 ^ 1 ' ^ere (Vi : qi e q? d ® & B vk) A (i i' ^ Qi =t Qi>)- (28) 

i i 

These superpositions can consist at most of N Bloch waves with corresponding N different q- values q\,...,qx G 
q^ d © ObvK- The number of different bands involved in this superposition can be any number between 1 and N . 

Superpositions of the form given in Eq. (p?5|) have a special property following directly from Eq. : The expectation 
value of a periodic differential operator restricted to a finite domain of N periods [cf. Eq. (|2"Tj) ] is equal to the sum 
of the weighted expectation values of that operator with respect to the individual Bloch waves constituting the 
superposition: 

(A) Nd = J2^K^) Nd = y E(<xiK,«\ A \°>iK,«)™ = E \<*\ 2 (K,«\ a \<,«)n* ( 29 ) 

i j i i 

Note that these expectation values are time independent. Indeed the time dependence of the wave function |\&(t)) = 
Oii(0) exp(— tE(m, qi)t/K)\^^. q .) is not reflected in the expectation value: 

(A(t)) Nd = ]T k(o)| 2 << = (A(o)) Nd . (30) 
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A. A set of operators which have real expectation values with respect to Bloch functions in C 2 ([0, Nd]) 

Next we concentrate on a set of operators Ah that are known to be self-adjoint with domains that are dense in 
C 2 (R). First we note that such operators are in general non-self-adjoint over C 2 ([0, Nd]). As an example consider the 
momentum operator Ah = P = — on the domain [0, Nd]. Partial integration gives 

l-Nd 
JO 



-ih 



The operator is clearly non-symmetric (and therefore non-self-adjoint [01]) as the boundary term on the right hand 
side does not vanish for arbitrary Bloch waves. Therefore such operators do not have real expectation values with 
respect to arbitrary linear combinations of Bloch functions. 

The set of operators which we study is built by (a) operators of the form Djj = (— i4-) n , n £ IN, (b) real valued 

periodic functions B : R — > R,B(x) = B(x + d), (c) symmetrized products of operators and B, e.g. Ah = 

DjpBDg\ Ah — BD^B, Ah = \ (d^B + BD^^j and similar products with more operators. Such operators 

appear frequently in the context of periodic semiconductor heterostructures (superlattices). For a periodic potential 
V(x) and a periodic effective mass m(x), the momentum operator p, the velocity operator v, and the most commonly 

used Hamiltonian H are given by p = hD$ , v = § (b^m^ix) +m- 1 (x)D$\ H = ^-D ( ^ ) m- 1 (x)D { ^ ) + V(x), 

respectively. 

We now prove that the expectation values of Ah with respect to Bloch functions in C 2 ([0, Nd]), i.e. 
{^^ q ]An \^nq)Ndi are still real numbers. To this end we start with the first order differential operator D^' — 
Inserting (x|\I/^. q .) = (x|\f^ q ) = (x\^^ q ) — u n , q (x)) exp(igx) into Eq. (f3"Tj) leads to a purely periodic function 
u* n q u n _ q (x) in the boundary term, which therefore vanishes. Thus 

(K^K,)** = (K^H^KXm (<. q \D H ] \< q )N d e R. (32) 

For an n-th order differential operator Ah iterative partial integration leads to n boundary terms. Regardless of 
the order n the exponential terms exp(jfiqx) cancel and purely periodic functions remain again. Thus 

(Ht* v \A H \** q ) Nd eR. (33) 

As a consequence, together with Eq. ([29]) , the operators which we studied in this section have real expectation values 
for superpositions of the form of Eq. (|2"5|) , i.e. 



(Vt : qi € q^ d © Qbvk) A (i ^ i' q, ^ q V ) => 



B. Resonant superpositions 

Choosing one of the qi's in the superposition Eq. (|28p to be a resonant Bloch wave number qr es , results in a linear 
combinations of resonant states composed of either even or odd k values only: 

m = J2 ak K k M } ' Eki 2 = 1 ' ( 35 ) 

kefi ' k 

where Si depends on whether N is even or odd: For N = 2m we have the set of possible even k values Si = {k 6 
2Z, -N + 2 < k < N} and the set of possible odd k values Si = {k e 2Z+1, -N + 1 < k < N - I}. For N = 2m + 1 
we have Si= {k e 2Z, -N + 1 < k < N - 1} and Si = {k E 2Z + 1, —N + 2<k< N}, resp. 
Then 

(A) Nd = E I«*ft** )flffl l4** e R - ( 36 ) 
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Note that we have implicitely included the band edge functions \^f B (0 > ), , m ) in the decomposition which are 

n k ,qi. e ' B ' n k ,q}. es " 

strictly speaking non-resonant, but their corresponding Bloch wave numbers belong to the extended group O ros - 
For a CQC system we have seen that (cf. Section lilBI) Qqqq = Sres and that \af-\ = |ck Consequently 

(A) Nd = Yl l«»=| 2 (K «=> <*> )^+(*f a^\ A K a™)**) E R > ^ 

^— ' V n k,q C QC n k,1cQC n k,q C QO n k,1cQC I 

l<k<N— 1 



and operators A that fulfill 



<tf B w w ) Nd = -(^ B (38) 



have necessarily a vanishing expectation value (A) 



Nd- 



C. Velocity Expectation Value at Resonance 



Equation (|36|) shows that the velocity expectation value for a linear combination of several resonant Bloch waves, 
Eq. (|3"5j). is given by 

(i)Nd = ^ \ a k\ 2v g( n k,q { r k el) & R, (39) 

where w 9 (n., fc) = hr^dEn ^j dq denotes the group velocity, which is the velocity expectation value of a Bloch wave. 
For a superposition of two energetically degenerated resonant solutions only, i.e. 

|*)=a + |*^ s )+a_|^_ g ( S ), |a + | 2 + M 2 = l, (40) 
Eq. (|39[) and f s (n, — g) = —v g (n, q) show that the velocity expectation value is simply given by 

(v) Nd = {\a + \ 2 -\ a ^)v g {n,q^ s ) 6 R. (41) 

This result has been obtained in Ref. [l(| by using rather elaborate calculations. For a CQC systems, where \a+\ = 
a- 1, (v)pfd — is obviously fulfilled as it should be. 



VI. CONCLUSIONS 



In this paper we have discussed different properties of Bloch waves in finite periodic systems with an integer 
number of unit cells related to their Bloch wave numbers. We have summarized wave number quantizations and 
their mutual interrelations for finite periodic systems in three different cases: periodic boundary conditions [PBC, 
Eq. iJSJ)], complete quantum confinement [CQC, Eq. ©J and resonant states in open systems [OBC, Eq. Ijlip], Further, 
useful vanishing integral identities [cf. Eqs. (|13|) - (|17p ] for open finite periodic systems, which are generalizations of 
well known similar identities for PBC systems have been derived. The main innovation in this respect is related to 
the fact that symmetry properties of the PBC system manifested by the special role of the set 0BvK [cf. Eq. ©] 
reappear in finite systems despite their lack of obvious symmetry. We have shown that all matrix elements of periodic 
differential operators [as defined in Eq. (|2ip ] - restricted to a finite domain with an integer number of unit cells - for 
two different Bloch waves (of the same or of different bands) with wave numbers in the same coset of ObvK vanish. 
Consequently these two Bloch waves must also be orthogonal. Finally we have considered superpositions of Bloch 
waves with wave numbers that belong to the same coset of 0b v k- The expectation value of the discussed periodic 
differential operators [cf. Eq. I|21[)] with respect to these superpositions is time independent, real-valued and equal to 
the sum of the weighted expectation values with respect to the individual Bloch waves constituting the superposition. 
In addition some specific applications of these general results are presented. 
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APPENDIX A: GROUP THEORETICAL PROOFS 

We start by introducing some appropriate general group theoretical concepts and then use the symmetry properties 
of the problem at hand for performing the proof. 

1. Relevant groups and irreducible representations thereof 

The symmetry group of (one-dimensional) infinite periodic systems with period d is the one-dimensional discrete 
translation group Td, which is generated by an elementary translation d: 

T d = dT, n, m G Z, g n = nd, g m = mdeT d ^ g n * g rn = (n + m)d = g n+m . (Al) 

As %i is abelian it is well known that its irreducible representations are one-dimensional and are of the form 

D T d i9n) = <ixp(iqnd), q e (— w/d, n/d]. (A2) 

The irreducible representations of Ti form an abelian group Ti with the tensor product as group operation: 

D T d (9n) D T d (9n) = expand) exp(iq 2 nd) = exp^^ + q 2 )nd] = exp^^ © q 2 )nd] = 

= D^( 9n ), (A3) 

where q\ © q 2 denotes the minimal residue of (q\ + q 2 ) mod 2n/d, or in physical terms, the sum qi + q 2 folded back 
into the first Brillouin zone (— n/d, ir/d] [cf. Eq. (0)]. 

Indeed, (i) the tensor product of two different irreducible representations of this group is also irreducible, (ii) this 
tensor product is associative and commutative, (iii) the unit element is given by the trivial representation = f , 

whereas (iv) the inverse element of T)\ is its complex conjugate . 

Clearly the group Td and [cf. Eq. (|12p ] are isomorphic through the mapping q <-* exp(iq o). 

We now introduce the symmetry group of a system having translation symmetry Ti in which additionally periodic 
BvK boundary conditions are imposed. This symmetry group is the cyclic factor group T per = Td/T^d, where T^d is 
a subgroup of Td generated by the elementary translation Nd. 

Tper is explicitly defined as 

Tper = {h = 0,h 1 = d,h 2 = 2d,..., h N -i = (N - l)d}, 
hiShj = (i + j)d mod Nd, (A4) 

where EH is the group operation in T per . All translations of the form mNd, m G Z act trivially within T per . 
The irreducible representations of the cyclic group T per are given by 



2kir 

D q T k (hi) = exp(iq k ld), q k = ^j, < k, I < N - 1. (A5) 

Nd 



Analogous to the case of the group Ti the irreducible representations of T per form an abelian group T per which is 
isomorphic to ObvK- As a side note we mention that - since T per is finite and abelian - it is also isomorphic to T per . 

The irreducible representations of %i with q = q-g^K are irreducible representations of T per . This follows from the 
fact that (i) the set of all possible indices q k of irreducible representations of T per coincides with the set QbvK of 

all <7 B y K , cf . Eq. (|6|) , and (ii) any general translation g n = nd G Td can be represented as nd — Id + mNd, where 
lid 'X'pgj- w 
Indeed 

Df; K {9n) = exp(iqW K nd) = exp(iq^ K ld) = D^fa). (A6) 
In summary, we have introduced a number of groups which are related as follows: 

{T Nd ,*} G (T d ,*), 

(T per , EH) = {T d , *)/ (T Nd , *> = (Tper, 6§> = (QbvK, 8) C (Q d c , ©) = (f d , (9)>. (A7) 



9 



2. "Selection rules" theorem 



Consider a function f(x) = J2k a k^k(x), which transforms according to an irreducible representation D of some 
group G, i.e. 

/Or 1 *) = E feJ D(g) kJ a^ k (x), Vg € G. (A8) 
It follows from standard representation theory that 

D^D° => [ dxf(x)=Q, (A9) 
Jm 

where D° is the identity representation of G and Ai is a group invariant manifold. This general theorem [l5|, [l6| is 
the basis of the quantum mechanic selection rules and follows from the well known orthogonality relations of group 
representations [ia. [ltl Il7| . 

In the case of T per all irreducible representations are one-dimensional, i.e. f(x) ~ exp(iq^ K x) and M = [0, Nd}. 
In this case one can directly see the validity of Eq. (|A9j) : 



rNd /N-l x - a 

fc + ^ ^/ per (x)e^BvK- =\Y J cxp(zq^ K ld)\ ^ / per (x)e^Bv K - dx 



1 — cxp(±i27rfc) 
1 -exp(±i2nk/N) J 



f per (x)e iq ^ x dx = 0, (AlO) 



due to the vanishing numerator. This holds only if k is not a multiple of N, which is, of course, fulfilled for k £ A'bvK 
[cf. Eq. ©]. 



3. Proof of Equations (|13H15p 

To prove Eq. (fT3)) one should take into account that provided a function f(x) transforms according to an irreducible 
representation of a translation group its derivatives 4-^f{x) also transform according to the same representation of 
this group. This fact follows from the invariance of the derivative operator with respect to translations. 

The integrand of Eq. (|13j) transforms with respect to the infinite translation group 7^ as 

D a rd ®\^D^=D a rd ®D®^ =D q ™, k^O. (All) 

(fc) (k) 

From Eq. (IA6|) we know that -D^ vK = . This result together with Eq. (|A9j) completes the proof. 

Taking into account that [D^]* = D T ^ 3 , Eq. (flj) follows directly from Eq. (|A11|) . Eq. ((T5J) follows also directly 
from Eq. (|A11|) taking additionally into account the translational invariance of the differentiation operator. 

It should be noted that the essential step in the proof can be seen as a embedding the problem at hand into another 
one with appropriate symmetries. The result is then a consequence of the symmetry properties of this second problem. 
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